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Abstract. The construction of 'PT-symmetric quantum electrodynamics is reviewed. 
In particular, the massless version of the theory in 1 + 1 dimensions (the Schwinger 
model) is solved. Difficulties with unitarity of the ^-matrix are discussed. 



1. PTQED 

Quantum electrodynamics (QED) is by far the most successful physical theory ever 
devised [1] . However, although its reach includes all of atomic physics, there are a myriad 
of phenomena that we do not understand. Thus it is essential to explore alternative 
theories, in the hope that we may be able to describe aspects of the world that are as 
yet not under our understanding. 

One very promising new approach to quantum theories are those included under 
the rubric of non-Hermitian theories, in particular theories in which invariance under 
the combined operation of space and time reflection PT replaces mathematical Dirac 
Hermiticity in order to guarantee unitarity of the theory. (For a recent review, see [2].) 
Little work, however, has been done on applying this idea to quantum field theory. This 
paper represents our continuing effort to develop a PT-symmetric version of quantum 
electrodynamics, in the hope that a consistent theory, with a unitary S-matrix, results 
that may eventually find physical applications in nature. 



1.1. Transformation properties 

At the first International Workshop on Pseudo-Hermitian Hamiltonians in Quantum 
Physics [3] a "PT-symmetric version of quantum electrodynamics was proposed. A 
non-Hermitian but PT-symmetric electrodynamics is based on the assumption of novel 
transformation properties of the electromagnetic fields under parity P transformations, 
that is, 

P : E -> E, B -> -B, A -> A, A -> -A , (1) 
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just the statement that the four-vector potential is assumed to transform as an axial 
vector. Under time reversal T, the transformations are assumed to be conventional, 

T : E -> E, B -> -B, A -> -A, A — > A . (2) 

Fermion fields are also assumed to transform conventionally. We use the metric 
<T = diag(-l, 1,1,1). 

1.2. Lagrangian and Hamiltonian 

The Lagrangian of the theory then possesses an imaginary coupling constant in order 
that it be invariant under the product of these two symmetries: 

I 1 I 771 1 

C = --F^F, U - -V, 7 Vt^ - y ihV + ^7Ve#^. (3) 

Here, because we are discarding Hermiticity as a physical requirement, it is most 
appropriate to use a "real" field formulation, where correspondingly the (antisymmetric, 
imaginary) charge matrix q = 02 appears. Furthermore, 7°7 M is symmetric and 7 is 
antisymmetric. In the radiation (Coulomb) gauge V ■ A = 0, the dynamical variables 
are A and ip, and the canonical momenta are 7Ta = — E T , ir^ = ^ip, where T denotes the 
transverse part, and so the relation between the Hamiltonian and Lagrangian densities 
are 

H = E 2 + E • VA° + - C. (4) 

Then, if integrate by parts and use V • E = ij°, we find that the corresponding 
Hamiltonian is 

H = j (dr) {i(£ 2 + B 2 ) + ^ [ 7 V (yV fc - ieqA k ^ + m 7 °] ^} . (5) 
We can also obtain this same Hamiltonian from the stress tensor 

= --V>7° {Yd" + Yd 11 ) ip + F^F V X - - (fA v + fA") + g" u C. (6) 

1.3. Current density 

The electric current appearing in both the Lagrangian and Hamiltonian densities, 

f = ^tV^, (7) 
transforms conventionally under both V and T: 

Tf(x,t)T= ^.y*,-t). (8b) 
This just reflects the normal transformation properties of the fermion fields. 



VT -Symmetric Quantum Electrodynamics — VTQED 3 
1.4- Equal-time commutation relations 

We are working in the Coulomb gauge, V • A = 0, so the nonzero canonical equal-time 
commutation relations are (the fermion index a includes both the Dirac and charge 
indices) 

{^«(x,*),V&(y,*)} = <W*(x-y), (9a) 



[4-(*),^(y)] = -i 



On - 



<f(x-y). (96) 



V 2 

V • A T = V • E T = 0. (10) 
1.5. The C operator 

As for quantum mechanical systems, and for scalar quantum field theory, we seek a C 
operator in the form 

C = e Q V, (11) 
where P is the parity operator. C must satisfy the properties 

C 2 = 1, (12a) 

[C,VT\=0, (126) 

[C,H] = 0. (12 c) 

From the first two equations we infer 

Q = -POP, (13) 

and because PT = TP, 

Q = -TQT. (14) 

The third equation (112c|) allows us to determine Q perturbatively. If we separate the 
interaction part of the Hamiltonian from the free part, 

H = H + efli, (15) 
and assume a perturbative expansion of Q: 

Q = eQ 1 + e 2 Q 2 + ... J (16) 
the first contribution to the Q operator is determined by 

[Q 1 , J ff ] = 2 J ffi. (17) 
The second correction commutes with the Hamiltonian, 

[Q 2 ,H ] = 0. (18) 
Thus we may take 

Q = eQ 1 + e 3 Q 3 + ..., (19) 

which illustrates a virtue of the exponential form. The 0(e) term was explicitly 
computed for four-dimensional QED in 2005 [I]. 
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2. (PTQED) 2 

However, the above perturbative construction of C fails for 2- dimensional PT-symmetric 
QED. We are here discussing the Schwinger model [3 [6], [7J [8] . In two dimensions, the 
only nonzero component of the field strength tensor is F 01 = E, and the Hamiltonian 
of the system is 

c r 1 1 TTl 

H = jdx y-E 2 - ij x A x - -^Vdiip + y^Vj • (20) 

As before, we choose the radiation gauge because it is most physical: 

V ■ A = d 1 A 1 = 0, (21) 

and then the Maxwell equation becomes 

d 1 E 1 = -dfA° = if, (22) 

which implies the following construction for the scalar potential 



A\x) = -- dy\x-y\f(y). (23) 

L J— oo 

2.1. Construction of E 

Without loss of generality, we can disregard A±, and then the electric field is 

i r°° 

E(x) = - dye(x-y)f(y), (24) 

Z J —oo 

with 

f 1, x > y, 

e( x ~y) = \ °> x = V' ( 25 ) 
[-1, x< y. 

Thus the electric field part of the Hamiltonian is 

J dx^E 2 = - ^J dxdydze(x-y)e(x-z)f(y)f(z) 

= - \LQ 2 + \J 'dydzf(y)\y- z\f(z), (26) 
where L is the infinite "length of space" and the total charge is 



Q = J dyf(y). (27) 
As this is a constant, we may disregard it. 

2.2. Form of Hamiltonian 

Thus we obtain the form found (for the conventional theory) years ago by Lowell Brown 

m- 

H = - A J dydzf(y)\y - z\f{z) - J dx {^tV^ - yV^V} ■ (28) 

This resembles 4 theory, and for the same reason, we cannot calculate the C operator 
perturbatively. 
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2.3. Functional integral 

It may be useful to rederive the Hamiltonian using the functional integral approach. 
The partition function is 

Z = [[di>][d(j)]e i f dtdxL , (29) 

where in the Coulomb gauge (with A = <p) 

] 1 1 

L = g (<9i0) 2 - ij'V - Vt^ - y^V (30) 
We integrate out the scalar potential by completing the square, 

where 

,-5l*-»l- (32) 

The functional integral on <fi is carried out over real values of <p(x). Then, the partition 
function can be written as 

\d^}e i I dtdx{7Tljj - n \ (33) 



where the momentum conjugate to if) is 



The result for H, given in (1281) . is reproduced. Because the sign of the quartic term in 
H is reversed, presumably we can no longer regard [dip] in fl33l) as over "real" values of 
if). Since Grassmann integration is a formal procedure, it is not immediately clear how 
to proceed. Henceforth, we will set the fermion mass m = 0, so we will refer to the 
Schwinger model proper. 

2.4- ETCR of currents 

It is easy to check that 

\j°(x,t),j°(y,t)} = 0. (35) 

However, it requires a point-splitting calculation [which does not modify ([35]) ] to verify 
that 

\j%x,t)j\y,t)] = -^5(x-y). (36) 

The key element in the latter is that the singular part of the 2-point fermion correlation 
function is given by the free Green's function: 

{M^Wvh'h) = \G a p{x - y), (37a) 

G(z) = ~ (376) 
This agrees with the massless fermion propagator found in (15960 below. 
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2. 5. Conservation of electric charge 
The electric current is exactly conserved: 

d f = - i [f,H] = -d lJ \ 

or 

d,f = 0. 

2.6. Axial-vector anomaly 

In 2-dimensions, the dual current is 

Now, using the above commutator between j° and j 1 , we find 
d *f = d j 1 = - i [j lj H] 



+ - 
i 



ii(^)>7 / d y dz f(y)\y - z \f( z ) 



4 

so from (1361) this can be rewritten as 

2 

d,*f(x) =-^J dydzdj(x -y)\y- z\j°(z) 
= - ^° = -E. 

7T 7T 

This is the two-dimensional version of the famous Schwinger-Adler-Bell-Jackiw 

2. 7. Schwinger mass generation 

Combine the current conservation and axial-current non-conservation: 

di[d j° + dd 1 = 0} 
ie 2 

do[d j 1 + d l j° = —E}, 

IT 

together with the Maxwell equation (1221) to obtain (<9 2 = — <9q + d\) 



7T 



This corresponds to a spacelike singularity, a pole at 



e = -d 2 = 



,2 
7T 



implying complex energies! 
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v, k 



—k — p 

Figure 1. Lowest-order vacuum polarization graph 
2.8. Perturbation theory 

This result is consistent with perturbation theory, where in general we expect all we 
have to do is replace 

e -> ie. (45) 

In fact, the Schwinger mass comes from one-loop vacuum polarization. In particular, the 
C operator appears to have no effect on the weak-coupling expansion [TO] . In general, 
it appears only ephemerally in the functional integral formulation [TTj . 

Let us sketch the Feynman diagrammatic argument for the mass generation 
mechanism in the modified Schwinger model. That is, we calculate the vacuum 
polarization operator for massless QED in lowest order as shown in figure [TJ This 
corresponds to 

n"" = "(ie) 2 tr / ^y=?7- ; ( TP> , (46) 



{27r) d 1 p 2 1 m 2 + {p-k) 2 ' 

where we have regulated the integral by working in d dimensions. The trace is evaluated 
as 

tr^W/ = d (g aP g 1& + £°V 7 - # a V 5 ) • (47) 
The denominators are combined according to 

I poo rl 

dss due~ isx , (48) 



p 2 (p — k) 2 Jo Jo 
with 

x = (1 - u )p 2 + u(p - k) 2 = (p - kuf + k 2 u(l - u). (49) 

Now the integration variable is shifted, p — > p + ku, and odd terms disappear upon 
symmetric integration, leaving us with 

d d p 



roo roo 

IF" = -(ie) 2 d / dss / du 



o Jo J (2n) d 
- cT (p 2 ~ k 2 u{l - u 



2p»p u - 1k^k v u{\ - u) 

(50) 

Note that in d = 2 dimensions, the contraction (trace) of the tensor, IP^, vanishes, 
which indicates that, apparently, the only gauge-invariant result could be zero. That 
this is incorrect is the result of the quantum anomaly, which appears only by setting 
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d — 2 at the end of the calculation. Proceeding onward, we use the following momentum 
integrals: 

' T 



J d d pe~ sp2 = (/dp. 
J d d pp^p v e~ 



-sp 



1 d 



-sp 



1 (Tl\ d / 2 



(51a) 
(516) 



When the s integrals are now carried out, we obtain an explicitly gauge- invariant form: 



IP" = -2 1 ^(ie) 2 7r- d / 2 rfr(2 - d/2) f 1 du \k 2 u(l - u) 

Jo 1 

which is divergent for d — 4, but finite for d — 2: 



d/2— 1 



IP" 



7T 



k 2 



(52) 



(53) 



We can similarly calculate the correction to the fermion propagator, given by the 
mass operator, using the covariant photon propagator, 



ie 



d d k 7 M 7(p — k)Y 



kfjkv 



{2ll) d k 2 {p-k) 2 

The gamma matrix structures are reduced as 

7 A 7(p - A;)7 A = (d - 2)7(2? - k), 

'jk'j(p — k) 7/c = k 2 ^(p — k) + 7fc[(p — A;) 2 + A; 2 — p 2 ], 
and so the integral, after the substitution k — > A; + pu, can be evaluated as 



(54) 

(55 a) 
(556) 



IP 



ds s du 



71 
S J 



\ d / 2 



(d-3)(l-u)-u-p 2 — 



-su(l — u)p 2 



d/2-2 



7,^(2 - ,/2)(, 2 )^- 2 Ll ( ^lfL /2 " ^ 



0, 



(56) 



Vr(d/2 - i/2)2 rf - 5 / 2 r(d-i), 

as d — > 2, that is, the mass operator vanishes. 

In the conventional theory, the vacuum polarization, iterated, yields the boson mass 
generated in the Schwinger model: 



D(k) 



1 e 2 1 



+ . . . 



1 



k 2 k 2 n k 2 k 2 + e 2 /n'' 

while there is no correction to the fermion propagator, 

\G{p) = —. 
IP 

The corresponding iterated one-loop propagators in coordinate space are 
D ( x ) = ^- K o( mx ), 



G(x) = 



1 7 M ^ 
2n x 2 ' 



(57) 

(58) 

(59 a) 
(596) 
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where in the conventional theory m 2 = e 2 /ir, while the mass-squared is reversed in 
sign in the VT theory. These propagators are given in terms of the Euclidean distance 
x = yx 2 . Of course, in other gauges, there are corrections to the fermion propagator. 
However, in Schwinger's words |6j, "the detailed physical interpretation of the Green's 
function is rather special and apart from our main purpose." 



3. Conclusions 



Perturbation theory evidently fails to give a positive spectrum to the massless VT- 
symmetric electrodynamics in 2 dimensions. More generally, this reflects the lack of 
unitarity of the S-matrix for the theory. This phenomenon has already been noticed 
by other authors: For example, Jones [12] shows for localized non-Hermitian potentials 
probability is not conserved, unless the Hilbert space metric is changed; and Smilga 
[T3] shows that the ^-matrix of a non-Hermitian theory defined in terms of "transition 
amplitudes between conventional asymptotic states is not unitary." 

Therefore, it seems that nonperturbative effects (strong field effects) presumably 
resolve this issue. It is necessary to do more than merely compute the Q operator in 
field theory, which as we have seen makes no explicit appearance in the functional or 
perturbative formulation of the theory, but we must determine the asymptotic states, 
presumably defined dynamically. This is a formidable problem. 

Previous work on these questions has concentrated on quantum-mechanical and 
scalar-field theory examples. Clearly there are issues unsolved relating to fermions and 
gauge theories in the PT-context. In particular, the formal technique of Grassmann 
integration needs to be redeveloped. One must find ways to reformulate usual field 
theoretic tricks, such as bosonization. 

Another illustration of the failure of perturbation theory to capture the essence of 
a VT theory is discussed in the Appendix. 



Appendix A. Zero-Dimensions 

We contrast the zero-dimensional partition functions for a conventional and a VT- 
symmetric x 2+N theory: 

Z C N (K)= [°° dxe- x2 ~ 9x2+N - Kx , (1.1a) 

Z N (K)= [ dxe- x2 ~ 9x2{ix)N - Kx . (1.1b) 
Jc 

Such examples were discussed earlier in [14|. The contour C in the latter integral is 

taken in the lower half plane, along Stokes wedges centered on the lines 

Ntt Nix , . 

6 = , 9 = 7i , 1.2 

4 + 2A' 4 + 2A' v ; 

which have width 2tt/(A + 2N), so that the integrand decays exponentially fast. 
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Note that the PT-symmetric theory has a perturbation series that doesn't appear 
to know about the path of integration: 

7V+2" 

/ CI } 

9 



Z N (K) 



it exp 



• d V 



K 2 /8 



n=0 



-i)M n 



21+^/2 



1 AiT 

n\ Dn{N+2) I 71, 



1.3) 



where D m (x) is the parabolic cylinder function. 

When we set N = 2 and K = we get the —x 4 theory without sources, for which 
we have the closed form for the vacuum amplitude 

(I \ , /I s 



Z 2 (0) 



71 



±V9 



-l/8g 



'■1/4 



*9, 



1-1/4 



*9, 



;i-4) 



The conventional theory in the same situation has the closed form 

Zm = JLttoK* (1) . (1.6) 

Directly from (jl.4p . or from the previous expansion fll.3j) . we find the weak-coupling 
expansion (g — > 0) 



zm - x " ( i + 1'/ + ^|// 2 + • • • 



(1.6) 



the expansion of Z| differs only in the sign of g. The conventional theory is Borel 
summable, while the VT series is not. 

However, the correspondence is not so simple in strong-coupling. Even the leading 
prefactors in the strong-coupling expansions are different: (g — > oo) 



zm 



Z 2 (0) 



71 



20 1 /4r(3/4) 



7T 



2c/ 1 /4r(3/4) 



i r(3/4) 

VsT(5/4) 

i r(3/4) 

4^r(5/4) 



+ . . . 
+ . . . 



;i.7a) 
(1.76) 
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